Electron-Phonon Coupling in Highly-Screened Graphene 
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Photoemission studies of graphene have resulted in a long-standing controversy concerning the 
strength of the experimental electron-phonon interaction in comparison with theoretical calculations. 
Using high-resolution angle-resolved photoemission spectroscopy we study graphene grown on a 
copper substrate, where the metallic screening of the substrate substantially reduces the electron- 
electron interaction, simplifying the comparison of the electron-phonon interaction between theory 
and experiment. By taking the nonlinear bare bandstructure into account, we are able to show that 
the strength of the electron-phonon interaction does indeed agree with theoretical calculations. In 
addition, we observe a significant bandgap at the Dirac point of graphene. 



The electron-electron and electron-phonon interactions 
are two of the fundamental interactions in many-body 
physics, giving rise to superconductivity, Mott-insulating 
behavior, and other collective phenomena. These phe- 
nomena are often studied in association with graphene 
not only because graphene is a simple system featuring 
two carbon atoms per unit cell [l[ but also due to the 
unique potential of this material. However, despite this 
apparent simplicity, there have been many difficulties in 
matching theoretical predictions to experimental studies 
of electron-phonon coupling in graphene The nature 
of these discrepancies is due to the way the electron- 
phonon coupling constant A is extracted from the ex- 
perimental data, and more specifically the way the bare 
velocity is determined. Within the Migdal-Eliashberg 
regime, electron-phonon coupling results in single phonon 
excitations that can be treated as perturbations to the 
bare band dispersion and leads to a renormalization of 
the group velocity with respect to the electronic bare 
band. The relative change of the renormalized velocity 
with respect to the bare velocity provides a measure of 
the electron-phonon coupling constant A. Angle-resolved 
photoemission spectroscopy has shown to be an invalu- 
able probe to extract this constant since it can directly 
measure the single particle spectral function and hence 
the renormalized velocity. However, the correct deter- 
mination of A rests on an accurate determination of the 
bare velocity, which is often done by assuming a linear 
band approximation between the Fermi energy and high 
energy, a procedure that has been found to be grossly 
inappropriate for graphene [3]. 

The LDA band velocity is in fact known to change sig- 
nificantly over the relevant energy scales, which greatly 
affects the measured values of the electron-phonon cou- 
pling constant A if a linear band is assumed. Further com- 
plicating the analysis, electronic correlations are known 
to renormalize the bare-band velocity in a nonlinear man- 
ner, to a degree determined by the dielectric screening of 
the substrate It should also be noted that the 

electron-phonon coupling strength may be enhanced by 



interplay between electron-electron and electron-phonon 
interactions [8l-[lo| . Since the experimental bare band is 
so difficult to determine, ARPES studies typically resort 
to the linear bare band approximation when determin- 
ing the electron-phonon coupling strength in graphene, 
resulting in an over- or under-estimate of the actual 
electron-phonon co uphng constant when extracted from 
the real self-energy Ill4l7l| . 

In light of these difficulties, one way to simplify the 
study of electron-phonon coupling in graphene might be 
to grow graphene on a metallic substrate, a growth tech- 
nique that has recently become popular due to its rel- 
evance for technological applications [3] ■ On a metallic 
substrate, the electron-electron interaction in graphene is 
expected to be highly screened, which would remove ve- 
locity renormalizations due to electronic correlations and 
cause the bare dispersion (experimental minus electron- 
phonon interaction) to converge to the LDA result [l9l - 
2l| . In this highly screened limit, the curvature of the 



graphene LDA band structure may be taken into account 
when analyzing electron-phonon coupling. Therefore the 
presence of a metallic substrate allows us to examine the 
electron-phonon interaction with an accuracy unmatched 
in other systems, leading to a straightforward analysis of 
the experimental data. 

Here we present a high-resolution ARPES study of 
graphene grown on a copper substrate. Starting from 
a basic characterization of this system, which has never 
been studied before by photoemission spectroscopy, we 
observe sharp dispersions due to the copper substrate 
and graphene overlayer, including a band gap at the 
Dirac point of graphene. Proceeding to examine the 
many-body physics in highly screened graphene, we find 
an overall agreement between the experimental band- 
structure and the LDA band calculations. Taking the 
curvature of the LDA band into account, we find close 
agreement between experimentally extracted electron- 
phonon coupling constants and theoretical calculations, 
providing the first real measurement of the electron- 
phonon coupling constant and providing closure to a 
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FIG. 1: (Color online) (a) Partial map of the Fermi surface. 
Arrows and circles correspond to bands from the copper sub- 
strate and graphene overlayer, respectively. Dashed black 
lines correspond to the hexagonal Brillouin zone of graphene. 
The K point is labelled, while the F-point is not shown, lo- 
cated at (kxjky) = (0,0). The black horizontal line through 
the K point illustrates the orientation of the data taken in 
panels (b) and (c), while the vertical purple line illustrates the 
orientation of the data in panel (d). (b) ARPES dispersion 
and (c) EDCs taken along the F-K direction at ky = A~^, 
showing that the graphene bands are n-doped with a bandgap 
and intensity minimum at the Dirac point. The presence of a 
bandgap creates two peaks in the EDCs at the K-point (peak 
positions marked in red), (d) EDCs taken through the Dirac 
point along constant k^,. In contrast to panel (c), where pho- 
toemission matrix elements suppress one branch of the cone, 
the photoemission intensity in panel (d) is symmetric and 
allows the presence of a bandgap to be easily seen, (e) Angle- 
integrated spectra of panel (b) and of highly-doped graphene. 



long-standing debate. 

Samples were grown on copper films as previously re- 
ported High-resolution ARPES data were taken at 
BLIO.0.1 and BL12.0.1 of the Advanced Light Source at a 
temperature of 15°K after annealing samples to 1000°K, 
using a photon energy of 50eV. The vacuum was better 
than 3 x 10^^^ Torr. Potassium was deposited in situ 
with an SAES potassium vapor source. 

Figure 1 (a) shows a Fermi surface map with bands due 
to the copper substrate and graphene overlayers. Two 
sets of copper bands and Dirac cones can be distinguished 
due to the presence of rotated crystallographic domains 
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FIG. 2: (Color online) (a-c) ARPES dispersions for several 
dopings, taken along the F-K direction. Arrows indicate the 
Dirac point in each image. Panels (a), (b), and (c) correspond 
to dispersions a, 7, e, respectively, in panel (d). (d) Exper- 
imental dispersions (red) for several dopings extracted from 
MDC peak positions, and LDA bands (blue) for the same 
doping along the F-K direction. Greek letters a, P, 7, 5, e 
label these dispersions in order of increasing doping. Inset: 
Comparison of dispersions a and e, showing that the electron- 
phonon kink is stronger for the more highly doped dispersion. 
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FIG. 3: (Color online) (a) Experimental coupling constants 
are given as a function of electronic charge density in red. 
The data agrees with theoretical calculations 30], shown in 
black. The error generated when the linear band approxima- 
tion is applied to the curved (bare) LDA band is given in blue, 
for comparison, (b) A zoomed-out version of panel (a), also 
showing results from the cited references in green (different 
references have different symbols depending on substrate). 
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of the substrate. The Dirac cones of graphene are visible, 
and are electron-doped due to their proximity to the cop- 
per substrate. Although doping can change from sample 
to sample, typical values are approximately 2 x lO^'^cm"^. 
ARPES data through a single Dirac cone is shown along 
the F-K direction in figure 1(b). In this measurement ge- 
ometry, the photoemission intensity is suppressed along 
half of the cone 22 1 . 

The dispersion in the vicinity of the Dirac point has 
been the subject of some controversy in the past. The va- 
lence and conduction bands are not coUinear, possessing a 
region of vertical intensity between them. Whether this is 
due to the presence of a bandgap in the bare disgersion or 
a many-body effect has been hotly debated 23|,|2j|. In the 
present case of graphene on a copper substrate, the di- 
electric screening of the highly conductive substrate rules 
out the possibility of electron-plasmon coupling 1^ 2^ , 
and instead implies the presence of a bandgap at the 
Dirac point [l^ 26 1. We find the separation between 



valence and conduction bands to be somewhat sample- 
dependent, with a typical bandgap of 400±50 meV, when 
determined from the separation between the peaks of en- 
ergy distribution curves (EDCs, intensity profiles at con- 
stant momentum), as shown in figures Ic and Id, with 
two peaks visible at the Dirac point momentum. The 
angle-integrated intensity (figure le) shows a V-shaped 
intensity profile, with a minimum at the Dirac point en- 
ergy, and increasing intensity away from the Dirac point. 
Far from the Dirac point, the valence and conduction 
bands are not coUinear, with an overall offset of 100±30 
meV, suggesting an unusual band gap opening mecha- 
nism [27|. This behavior is clearly similar to results re- 
ported previously, although the size of the bandgap at 
the Dirac point in this sample is larger 23, 2^. In figure 
le, a dip can also be seen at the phonon energy, where 
the density of states is renormalized due to electron- 
phonon coupling 28, [i^. In the past, calculations of 



electron-p honon c oupling have differed from experimen- 
tal results 
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15l | . This disagreement is believed to have 
three sources: first, the bare band is not completely lin- 
ear and may have a positive or negative second deriva- 
tive depending on the magnitude of the screening of the 
electron-electron interaction 0-0] and the direction along 
the brillouin zone second, the electron-electron inter- 
action is believed to enhance the electron-phonon cou- 
pling strength [il; third, the finite resolution of the ex- 
periment may lead to some error in the extracted band 
velocity [l^l . The first two effects might be eliminated by 
the presence of a metallic substrate as this is expected to 
screen the electron-electron and electron-plasmon inter- 
actions, and also to eliminate the enhancement effect. In 
the limit of infinite screening, metallic substrates are ex- 
pected to cause the graphene dispersions to converge to 
LDA calculations. This greatly simplifies the extraction 
of the electron-phonon self-energy. 

The electron-phonon interaction is visible in the pho- 



toemission spectrum in two ways: the real part of the 
self energy ReE modifies the band position; while the 
spectral width of the bands is proportional to the imagi- 
nary part of the self energy ImE [3l| . Peak positions and 
peak widths can be obtained by fitting lorentzian peaks 
to the momentum distribution curves (MDCs), the inten- 
sity at constant energy as a function of momentum. The 
strength of the interaction, given by the coupling con- 
stant A, can be extracted from either part of the self en- 
ergy, although in practice the real part of the self-energy 
is often more reliable, since the imaginary part is more 
sensitive to noise and the influence of impurity broaden- 
ing. We have therefore focused on the real self-energy in 
our analysis. 

Knowing the bare graphene band, ReS is given as the 
difference between the experimental and bare band po- 
sitions. From knowledge of ReE, the electron-phonon 
coupling constant A can be expressed as 



At — — 



dReY.k{E) 



or equivalently. 



A* 



dE 



Vk{Ef 



E — Ep 



(1) 



(2) 



where v'^{Ep) and Vk{Ep) are the bare and renormalized 
velocities at the Fermi level, respectively. However, the 
bare band of graphene is not linear, so the method of 
extracting A according to the formula 



V2 



1, 



(3) 



(where vi and V2 are the band velocities at higher and 
lower binding energy than the phonon, respectively) does 
not work, nor will any other method that assumes a linear 
bare band. (21 

The extracted coupling constants are compared with 
electron-phonon coupling calculations [s^l in figure 3. 
The agreement between experiment and theory is strik- 
ing, providing the first experimental support of theoret- 
ical electron-phonon calculations. Having said this, our 
analysis may require a small correction. Ab initio calcu- 
lations expect a finite el-ph self-energy even at high bind- 
ing energies 32|. This differs from our results, where the 
LDA band gives good agreement with experiment at high 
binding energies. This discrepancy could correspond to a 
difference in the measured coupling constant of approx- 
imately 0.015±0.005, and may derive from two physical 
origins: 1) It is likely that the metallic substrate does not 
perfectly screen the electron-electron interaction in the 
graphene overlayer. Since the electron-electron interac- 
tion increases the band velocity and the electron-phonon 
interaction decreases it, it is possible that at high binding 
energy both renormalizations affect the band velocity by 
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similar amounts and essentially cancel, leaving the exper- 
imental velocity to agree with LDA. 2) It is also possible 
that the LDA band is not a perfect description of the 
bare band dispersion, due to the presence of the gap at 
the Dirac point. 

It should also be noted that electron-phonon coupling 
has been studied on a metallic substrate in the past. 
One study of graphene on iridium attempted to extract 
the electron-phonon coupling constant in a self-consistent 
manner, but obtained a surprisingly lar ge v alue due to 
the approximation of a linear bare band [13| . 

To illustrate how much of a difference the linear bare 
band approximation makes, we have also applied the lin- 
ear approximation to just the curved LDA band (which 
does not include electron-phonon coupling), where we 
take to be the slope of the line that intersects the 
LDA band at energies E = Ep and E = Ep-^A eV using 
Eq. [2] The results, shown as the blue "LDA" line in fig- 
ure 3, correspond to the linear bare band approximation 
when no electron-phonon coupling is taking place. Data 
from previous works are given in panel (b). The error 
from the bare band approximation is more than twice as 
large as the actual electron-phonon coupling constant 2] . 
On the other hand, in cases where the bare band curves 
in the opposite direction ( "concave- up" ) , such as in the 
presence of strong electron-electron interactions [sI-ItI [33| . 
for the LDA band on the opposite side of the Dirac 
cone (along the M-K-Gamma direction) [J] , or for bilayer 
graphene [17|, the linear approximation underestimates 
the coupling constant. 

In conclusion, we have shown for the first time that 
the magnitude of electron-phonon coupling in graphene 
agrees with theoretical calculations. These results settle 
a long-standing controversy in the field, confirming the 
validity of theoretical calculations, and casting doubt on 
the conclusions of many experimental works. This work 
is also generally applicable to future experiments that 
require studying electron-phonon coupling in materials 
with nonlinear bare bands. We have also shown that 
there is a significant bandgap in graphene grown epitax- 
ially on a copper substrate, a discovery which may pave 
the way for future technological applications. 
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